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ABSTRACT 

We study systems of M Potts models coupled by their local energy density. 
Each model is taken to have a distinct number of states, and the permutational 
symmetry Sm present in the case of identical coupled models is thus broken initially. 
The duality transformations within the space of 2^^ - 1 multi -energy couplings are 
shown to have a particularly simple form. The selfdual manifold has dimension 
Dm = 2*^-1 - 1. Specialising to the case M = 3, we identify a unique non- 
trivial critical point in the three-dimensional selfdual space. We compare its critical 
exponents as computed from the perturbative renormalisation group with numerical 
transfer matrix results. Our main objective is to provide evidence that at the 
critical point of three different coupled models the symmetry is restored. 
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1 Introduction 



In the study of coupled models, and also of disordered models in their rephca formulation, 
the permutation group symmetry Sm is supposed to play an essential role ||, 
Namely, the interaction part of the action for a set of M identical coupled models 

Ant oc / d'^x gJ2^a{x)eb{x) (1.1) 

is explicitly invariant with respect to any permutation of the models. Here a and b are 
replica indices, and {ea{x)} designates the set of local energy operators. In the lattice 
definition of such coupled models the interaction part of the Hamiltonian takes a similar 
form, with only / d^x replaced by a summation and {^^(x)} by an appropriate lattice 
expression^. 

When one introduces asymmetric couplings, by generalising the common coupling 
constant 5^ to a matrix gab, 

Aint oc d'^x ^gabea{x)eb{x) (1.2) 

a perturbative renormalisation group (RG) analysis reveals that the Sm symmetry is 
restored at the fixed point. A detailed study of this scenario, within the e-type pertur- 
bative RG for coupled Potts models, was carried out in ||^. Supposing all components 
of gab to stay of order e, their initial values being all positive^], it was shown in |^ that 
the only non-trivial fixed point, having one attractive direction, all other directions being 
repulsive, is that with gab = This type of restoration of the symmetry Sm could be 
^In the case of the Potts model with nearest-neighbour spins (Tj-"'' and crj"'' in the rephca a, one has 

^In 1^ the interaction part of the action Ai^t (1^) was actuahy defined with an extra minus sign, so 

that initiaUy aU the components of gat were taken to be negative. 

•^Ref. also identified other fixed points, which are related to the basic one (with all gab equal and 

of the same sign) by changing the sign of some of its components. This is equivalent to switching the 

sign of a certain number of energy operators. This last operation is a symmetry of the individual Potts 

models, corresponding to their self-duality (note however such duality transformations on the individual 

models are not directly related to the global duality transformations on the entire coupled system to be 

discussed in section This argument implies that the critical properties of the various critical points 

classified in B should be equivalent. 
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called "soft universality". 

In this paper we are going to argue for a "strong universality" in the criticality of 
coupled models. We shall mainly be interested in coupling M = 3 different models^, via 
(|1.2|), with {Sa} belonging to Potts models with different number of states 5'2, ^s}- 
This breaks the permutational symmetry in a "strong" sense. Still, the RG calculations 
show the existence of a single fixed point with all {gab} being positive, like it is the case 
for identical models. 

To our knowledge, Simon |^ was the first to apply the RG analysis to a set of different 
coupled Potts models. The most general model studied by this author was that of Mi 
Potts models with qi states and M2 Potts models with q2 states (gi 7^ q'2), all of them 
being coupled. After determining the fixed point structure, he computed the dimensions 
of the spin operators, as well as the RG equations for the energy operators to two-loop 
order. The effect of disorder on these coupled systems was also analysed. 

Here we generalise the RG calculations of |^ to the case of three different coupled 
Potts models (gi 7^ ^2 7^ Qs)- We shall compute the dimensions of energy operators, with a 
special focus on the symmetry of the theory, at the non-trivial fixed point that generalises 
the one found in [g] for three identical models. Within the space of couplings {gab}, the 
new fixed point is stable in one direction and unstable in the others, the topology of the 
RG flows being similar to those of [0]. But there is also one apparent difference: The 
permutational symmetry has disappeared, the coupled models being different. 

The purpose of this paper will be to provide evidence, by using various methods, that 
at the fixed point of three different coupled models the apparently lost symmetry 5*3 is 
restored, implying a "strong universality" . 

This symmetry restoration cannot be observed on the level of the initial action ( |1.2|) , as 
discussed above, nor is it visible in the perturbative RG treatment, or in the Hamiltonian 
of the explicit lattice realisation. This is because are the energy operators of different 
models, with different scaling dimensions in particular. 

"^The duality transformations of section H are however valid for general M. 
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The way we shall check for the restoration of the symmetry is by looking at the spec- 
trum of scaling dimensions at the new fixed point, within the sector of energy operators. 
Like in the case of identical models [|I|, ^ ^ , the RG analysis implies that the three energy 
operators of the decoupled models {ei{x), e2{x), e^^x)} will rearrange as three particular 
linear combinations so as to form the new primary operators at the fixed point of the 
coupled models. 

In the case of identical models the corresponding linear combinations are easy to guess 
on symmetry grounds. The irreducible representations (irreps) of the group 5*3 in the 
basis {ei{x),e2{x),es{x)} consist of a (symmetric) singlet 

esix) = ei{x) + e2ix) + Esix) (1.3) 

and an (antisymmetric) doublet 

£:Ai = £i{x) - e2{x) 

(1.4) 

= - e-s^x) 

that act as the new primary operators at the fixed point ^ ^. The fact that the 
operators and belong to the same two-dimensional irrep means that their dimen- 
sions must coincide: A(£:ai) = ^{^A2)- These are however in general different from the 
dimension A(£:s) of the one- dimensional irrep. 

When coupling different models, the corresponding linear combinations of 61,82,8^ 
will have more complicated coefficients which have to be calculated by the RG technique. 
One will find something of the form: 

{el = KiiSi + Ki2e2 + i^ia^s 
el = K2iei + K22e2 + K2^ez (1.5) 

Since the initial dimensions A{ei), A{e2), A^e^) of the decoupled models differ, one may 
expect that the critical dimensions (RG eigenvalues) of the newly formed primary oper- 
ators el,€l,€'^ might all be different. 

Our argument is that, in the case of coupling different models, it is not the linear 
combinations ( |1.5|) which have to be examined to analyse the symmetry, but rather the 
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spectrum of their critical dimensions. Permuting 81,82-, 8^ in the combinations ( [1.5|) does 
not make much sense because they are different. To permute e]^,£2?^35 one first has to 
know their properties, their scahng dimensions, to decide if it makes sense or not. 

The conclusion will be that one has to study the spectrum of dimensions at the new 
fixed point. This provides a representation independent information, independent of the 
way one has defined the theory initially, as by its action (Hamiltonian) in Eq. ( |1.2|) . 

If the symmetry is restored then the dimensions 

A(5D, Hel). Hel) (1.6) 

should form a singlet and a doublet, as is the case when one couples initially identical 
models. 

The rest of the paper is organised as follows. 

In section |^ we analyse the case M = 3 by perturbative RG calculations. As will 
be explained towards the end of that section, in this case the RG results alone are not 
sufficient to decide whether the 5*3 symmetry is restored or not. We therefore propose 
to study the coupled models through a particular lattice realisation, following 0. To 
that end, in section ^, we extend the duality transformations derived in for symme- 
trically coupled models to the asymmetric case. The resulting relations are valid for 
any number M of coupled models, and for the most general asymmetric multi-energy 
couplings. In particular we determine the selfdual solutions, thus simplifying dramatically 
the subsequent numerical simulations. 

After recalling briefly, in section ^, the most efficient transfer matrix algorithm con- 
structed in p[, we show how it may be generalised to the case of asymmetrically coupled 
models. We then turn to the numerical analysis in section ^ We locate the non-trivial 
critical point on the selfdual manifold for various systems of M = 3 coupled models, and 
we provide accurate values of the central charge and the energetic scaling dimensions. 
We shall find strong evidence for a singlet/doublet spectrum at the new fixed point. This 
we interpret as a signal of the restoration of the 5*3 symmetry. 

Finally, section ^ is devoted to remarks and conclusions. 
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2 Renormalization group analysis 



In the continuum limit, the three coupled Potts models can be represented by the action 

A = j^Ai^^+A,^, (2.1) 

a=l 

3 

Ant = / d^xY.g',,ea{x)eh{x) (2.2) 

Here {^o'^''} represents three decoupled models, or more precisely the three decoupled 
conformal field theories for these models at their respective critical points. The interac- 
tion term Ai^t makes them coupled. In general, the initial couphngs {gab} are taken to 
be different. Since have different dimensions when the models are different, even 
if we start from identical couplings g^f, = go, they will become different in the course of 
renormalisation. 

We shall parametrise the dimensions of the energy operators {Sa} as in the papers 
i, 1, 1,0: 

(A(«))phy. = 1 - (2.3) 

where the physical dimension (A^"))phys corresponds to twice the conformal dimension, 
as usual. The quantity ea, which measures the deviation from the Ising model, appears 
also in the Coulomb gas parameter 

ic^^:^r = l-ea (2.4) 

For the Ising model one has e = 0, whence = | and (A£)phys = 1- The parameter 
is useful in analytic calculations because it is simpler to use than the central charge of 
the corresponding conformal theory. 

The details of the RG calculations are the same as in Refs. |^, 0. In particular, all 
the necessary integrals have been calculated in these papers. The generalisation to the 
case of three different models is straightforward, and we shall therefore only give the final 
results. 

It turns out to be convenient to redefine the expansion parameters {ea} and the 
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coupling constants {gab} as follows: 

3 



9a 



1 



-9ab 



(2.5) 



2 47r' 
From now on we shall adopt this convention, omitting the tildes throughout for simplicity 
of notation. 



The /^-functions for the interaction ( |2.2| ) are found in the form (to two-loop order): 

^ = Puiigab}) = eugu - ^i3te - Iguigh + gh) 

^ = Pniigab}) = ei3^13 - ^12^32 " |^13(^?2 + ^32) 
^ = /323i{gab}) = £23^23 " ^21^31 - |^23(^?2 + ^3l) 

Here ^ is the RG parameter, gab = gba by symmetry of the action 
defined eab = ea + eb. 

The renormalisation of the energy operators, to second order, is found to be given by 
the equations 



(2.6) 

, and we have 



In matrix form this reads 



b^a d^a 



d£a 

Kb 



;i - (-a)5ab - lab{.{gcd}) 



(2.7) 

(2.8) 
(2.9) 



/ 



(2.10) 



\{9l2 + 9h) -9i2 -gi3 

lab = -gi2 -U92I+923) -923 

^ -gi3 -g23 -^{gli + gh) ) 

To define the new primary operators and their scaling dimensions at the new fixed point 
we shall have to diagonalise the matrix Aab- It is however more convenient to regroup 
the terms —eaSab — lab so that 



A, 



ab 



Sab 
I 



A 



ab 



Kb 

^i-U9u + 9h) 

-gi2 
-gi3 



-gi2 -gi3 
^2 - ^{gli + gis) -923 
-g23 ^3-ligli + gh) 



(2.11) 



(2.12) 



and we need only diagonalise the matrix Kab- 

The non-trivial zeros of the /3-functions ( |2.(j| ) are found as 

9*12 = 7^13^23 {1 - |ei2 - ^(ei3 + £23)} 

9h = V^I^Il - |ei3 - i(ei2 + £23)} (2.13) 

^23 = Vei2ei3{l - |e23 - i(^12 + ^13)} 

They correspond to the non-trivial fixed point of the coupled models which we are inter- 
ested in. It is readily checked that the fixed point (|2.13| ) is stable in one direction and 
unstable in the two others, as is the case when one couples identical models. Furthermore, 
the topology of the RG flows is similar. 

Substituting the fixed point values of the couplings in the matrix Aab and diagonalising 
it, one obtains, after some algebra, the following expressions for the eigenvalues: 

a + b + c 3abc ^ ^. 

Ai = \ (2.14) 

2 a+h+c ^ ' 



A 



2,3 



a + b + c Qabc + + ab'^ + a^c + ac^ + + bc^ 
2 2(a + b + c) 



± — \ (6a^6^c^ - 2abc{a% + ab'^ + a^c + ac^ + + bc^ + + b^ + c^) 

2{a + b + c)^ ^ ' 

+ a%\a + bf + a\\a + cf + b\\b + cff' (2.15) 

where we have simplified the notation by means of the abbreviations a = ei2, 6 = £13 and 
c = €23. According to the definition ( |2.11| ) of Kab-, the dimensions of the new primary 
operators {^q}, cf. Eq. (|1.5|) , are related to the above eigenvalues through 

A(e:) = 1 - A, (2.16) 

for a = 1, 2, 3. We recall that all these calculations have been done to second order in e. 

It should be remarked that in the special case of two of the parameters being equal, 
b = c and a ^b, the expressions for the dimensions simplify considerably. One finds: 

AK) ^ (2.17) 

2 a + 2o 

Arf,^ = 1 — + ^ 2.18 

^ ^' 2 a + 26 a + 26 ^ ' 



At first order in e, the RG results for tlie dimensions 



. / *N a + 6 + c , , 

Aiel) = 1 + (2.19) 

A{e;) = A{e;) = i _ ^±|±^ (2.20) 

form a singlet and a doublet, in accordance with the scenario for the restoration of the 
5*3 symmetry discussed in the Introduction. The degeneracy of the doublet is however 
lifted at order e^, by the last term in ( |2.15|) . This is true even when two of the models 
are identical, and only for a = 6 = c does one recover the degeneracy ^{el) = A{el) 0. 

In section ^, where the results of our numerical work are presented, we shall show that 
this last result of the perturbative RG is wrongQ. We provide evidence that the splitting 
between the dimensions of eg and £3 is actually zero, and that the symmetry restoration 
scenario therefore holds true. 

We could suggest the following argument for the failure of the perturbative RG. 
The e-expansion calculations are valid for perturbed conformal theories because they 
respect the conformal symmetry, just as dimensional regular isat ion is valid in the context 
of pertubative calculations in a gauge theory because the method respects the gauge 
symmetry (otherwise the results would be dependent on the regularisation technique). 
In the present problem, the e-regularisation should be correct as far as the conformal 
symmetry alone is concerned. But in case of extra symmetries, such as S3, the method 
might well give wrong results. 

The RG formulae (|2.14|) -( ^.16 ) for the scaling dimensions, and for the central charge 



which will be given below, are still quite useful. Apart from the degeneracy issue just 
discussed, they compare well with the numerical results of section ^ when are small 
enough. 

As far as the dimensions ^{62), ^(£3), defined by Eqs. ( |2.15| )-( p.l6| ), are concerned, it 
is their mean value which compares well with numerical results. It is worth noticing that 
the splitting, namely the third term in (|2.15|) , is numerically smaller than the principal 

^We would still claim, of course, that our RG calculations are technically correct! 
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term, the second term in ( p.l5|) . This is because of the compensation of negative and 
positive terms in the third term of ( p^.l5D f| 



The central charge can be obtained in a simple way using Zamolodchikov's c-theorem 
P). This theorem provides us with a function of the couplings c{gi2, gis, g23) which 
decreases along the renormalisation flow and takes a value c{gl2, 9*3, 923) ^it the fixed 
point of the flow which equals the central charge of the associated conformal field theory. 

With the conventions ( |2.5| ) taken into account, the c-function is uniquely determined 

by 

dc{9i2, 913, 923) 3 

7^ = -7;Pab{9l2,9l3,923) (2.21) 

09ab 2 

c(0,0,0) = Co = Ci + C2 + C3 (2.22) 

where Cq is the total central charge of three decoupled models. 

From Eq. ( |2.6|) and Eqs. ( p.21|) -( ^.22| ) the c-function turns out to be given by: 



c(^i2,^i3,^23) = Co-^{a9l2 + b9l3 + C923) + ^9i29i3923 + ^{9u9l3 + 9u923 + 9I39I3) (2-23) 



At the fixed point we insert Eq. ( p.l3| ) into Eq. ( |2.23| ); up to the order the correction 
Ac of the central charge is: 

3 3 

Ac= — abc+ -(abc^ + bca'^ + cab'^) (2.24) 

4 8 



3 Duality transformations 

The possibility of endowing the Potts model with a duality transformation was one of 
the main motivations for introducing it [^]. The study of duality for several Potts models 
coupled by their local energy density was initiated by Domany and Riedel, who worked 
out the case of two models with qi and q2 states [|T0|. Technically, these authors used 
a particular version of the method of lattice Fourier transforms |jll||, due to Wu and 



^It is easy to check that the expression under the square-root sign of (2.15) is always non-negative, 
so that A2.3 are weh defined. 
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Wang [12|. Dotsenko et al. generalised the computations to three symmetrically coupled 
g-state models |0. However, it became clear that the complexity of the Fourier method 
grew rapidly as the number of models to be coupled was increased. Recently, it was 
pointed out by Jacobsen that trading the Potts spin variables for a formulation in terms 
of random clusters |TB[, or loops [Q, the duality relations simplified dramatically. This 



observation made it possible to work out the case of M coupled g-state Potts models, 
with the most general coupling by local energy densities consistent with an Sm symmetry 

i- 

We now show how to treat the even more general case, where each model does not 
necessarily have the same number of states, and the coupling by local energy densities is 
the most general one. 



3.1 General case 

Consider a set of M identical planar lattices £, which we imagine to be stacked on 
top of one another. On each lattice site i G £, and for each layer /i = 1,2,..., M, we 
define a Potts spin a-'^'' taking the values 0",^'^^ = 1,2,3, .. . ,g^. The layers interact by 
means of the reduced Hamiltonian 

'H = Y.'H^,, (3.1) 

where {ij) denotes the set of lattice edges, and the nearest-neighbour interaction is defined 
as 

n., = -Y.K,\{5(a\^\af). (3.2) 

i<z£ net 

By definition, the Kronecker delta function S{x,y) = 1 if x = y, and zero otherwise. We 
have here defined £ as the product set Y\'I^Li{^, C^^^}, so any one of the 2^^ subsets m 
S can be interpreted as a certain subset of layer indices. One can also think of i G S 
as specifying the state of a given edge (ij), meaning that in the interaction term ( p. 21) , i 
determines which of the layers contribute to the product Il^e^ of ^^e corresponding delta 
functions. The layers specified by some i G S then interact by means of the product of 
their local energy densities, through a coupling constant K^. 
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For later convenience, we shall represent a subset i G S as a list of M open (o) or 
filled (•) circles, the fith circle indicating respectively the absence (or presense) of the 
factor 6{a^^\aj^^). Furthermore, we sometimes use a single open (resp. filled) circle as 
an abbreviation of a list of M open (resp. filled) circles. 

When M = 1, the model defined by Eq. (^.21) reduces to the conventional Potts model. 



whilst for M = 2 it is identical to the Ashkin- Teller like model considered in Ref. [|10 
For M = 3, it is the asymmetric version of the model discussed in Ref. 0, which had 
q = qi = q2 = q3, and whose couplings possessed an 5*3 symmetry upon permutation of 
the layers. Using the above symbolic notation, this symmetry can be expressed by the 
identities 

K' — K — K — K 

K2 = Ko.. = K.O. = K,.o (3.3) 
= K.„. 



By means of a generalised Kasteleyn-Fortuin transformation the local Boltzmann 
weights can be recast as 



exp(-Hij) = n 



ic£ 



l + (e^^-l)n5(a(^\a 



(m) 
j 



(3.4) 



where we have simply used that exp{K6) = 1 + [exp{K) — 1]6, if 6 can only take the 
values and 1. Since furthermore 5^ = 6, expanding the product over i will lead to the 
equivalent form 



defining the coefficients hi. The normalisation of Eq. (|3^) is expressed by the fact that 
bo = 1. 



To relate the hi to the physical coupling constants Ki, we evaluate Eqs. ( ^.41 ) and 
( p.5| ) in the situation where S{al'^\ cxY^) = 1 for G and zero otherwise. This resulting 
equations 

exp(E^p|=E&P (3.6) 

\pci' I pee 
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can readily be solved by applying the principle of inclusion-exclusion [1^, yielding 



&. = E(-i)'''-'"'-P E^^P • (3-7) 

The partition function in the spin representation 

z = Y.E^M-n.i) (3.8) 

can now be transformed into the random cluster representation as follows. First, insert 
Eq. ( p.5| ) on the right-hand side of the above equation, and imagine expanding the product 
over the lattice edges (ij). To each term in the resulting sum we associate an edge 
colouring Q of where an edge (ij) in layer /i is considered to be coloured (occupied) 
if the term contains the factor a^'^), and uncoloured (empty) if it does not. In 

other words, an edge colouring Q is determined by specifying an edge state iij C S for 
every edge (ij) G C. 

The summation over the spin variables {a} is now trivially performed, yielding a 
factor of for each connected component (cluster) in layer fi of the colouring graph. 
Keeping track of the prefactors multiplying the (5-functions, using Eq. ( p.5|) , we conclude 
that 

^=En&f (3-9) 

where is the number of clusters in layer and is the number of edges (ij) G 
having the state i. 

It is worth noticing that the random cluster description of the model has the advantage 
that the only enter as parameters. By analytic continuation one can thus give meaning 
to a non-integer number of states. The price to be paid is that the are, a priori, non- 
local quantities. 

In terms of the edge variables the duality transformation of the partition function 
is easily worked out. For simplicity we shall assume that the couplings constants K£ are 
identical between all nearest-neighbour pairs of spins. The generalisation to an arbitrary 
inhomogeneous distribution of couplings is trivial; it suffices to let Ke depend on (ij) in 
Eq. (0). 
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By definition, we take the colouring state i C S of the edge (ij) G £ to be dual to 
the complementary colouring i* of its intersecting dual edge (ij) G C. In the symbolic 
notation introduced above, the complementarity operation * simply means replacing 
every • by a o, and vice versa. Also, note that we refer to dual quantities by a tilde 
throughout. 

To establish the duality transformations, we begin by postulating that the configu- 
ration Q, with all lattice edges coloured, be dual to the configuration Q* = Qo with no 
coloured (dual) edge. 

This requirement fixes the constant entering the duality transformation. Indeed, from 
Eq. (|3.9|) , we find that Q, has weight &f H^li Qfi, where E is the total number of lattice 
edges, and Qo is weighted by 6f Il^li <lf^, where F is the number of faces, including the 
exterior one. We thus seek for a duality transformation of the form 



where for any configuration Q the edge weights must transform so as to keep the same 
relative weight between Q and Q, as between Q* and Qo- 

An arbitrary colouring configuration Q entering Eq. ( |3.9| ) can be generated by applying 
a finite number of changes to Q,, in which an edge of weight b, is changed into an edge of 
weight be for some £ G S. By such a change, in general, a pivotal bond is removed from 
the colouring graph in each of some subset i' C i* of layers, thus creating new clusters 
in the corresponding layers. The weight relative to that of Q, will therefore change by 
(be/b,) Yl^^i' Qfj,. On the other hand, in the dual configuration Q a cluster will be lost in 
each of the layers {£')* fl £*, since each of the new clusters mentioned above will be 
accompanied by the formation of a loop in Q. The weight change relative to Qo therefore 
amounts to (b^* /bo) n/ie{€')*n^* Comparing these two changes we see that the factor 
n^e£' q^i cancels nicely, and the duality transformation takes the simple form 




(3.10) 



bi = {bi,/b,) W q^, 



(3.11) 



the relation with i = o being trivial. 



14 



Note in particular that Eq. ( p.llD with i = • imphes that b,b, = Y[^=i Qfi, since bo = I 
by definition. Then, duahsing Eq. ( |3.11D once again yields 



s=in?.-rn?.=&^ (3-12) 

so that the duality transformation is indeed involutive, as required. 

The duality relations ( |3.11| ) can be recast in an even simpler form by trading the 
random clusters for the loops surrounding them (and their duals) on the lattice 



medial to C Using the Euler relation, we find that Eq. (|3.9| ) must be replaced by 

z={e%] En^fn?^^ (3.13) 

\At=i / g ic£ fiee 
where, by a slight abuse of notation, we use the same notation Q for the loop and the 
cluster configurations, since they are in bijective correspondence. are now the number 
of closed loops in layer /i, and N is the total number of vertices in C. Note that the bond 
weights bi have now been replaced by 

Xi = bell<l>^'^'- (3-14) 



It is easily verified that the duality relations (|3.11 ) now simply read 



xe = Xi*. (3.15) 

This is our main result. 

In the case of a lattice C = C which is unchanged by the duality, such as the infinite 
square lattice, we can now search for selfdual solutions. These are obtained by imposing 
be = be, and read explicitly 

a;f<^-)=a;?/-). (3.16) 
Since Xo = 1 by normalisation, the selfdual manifold has dimension Dm = 2^^^^ — 1. 
Two special points always belong to the selfdual manifold. The first one is 

xe = 1, for any i C S. (3-17) 
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It is straightforward to verify that in terms of the original couphngs this means -ft'o - o.o - o = 
1 + (with the unique • at position fi) ioi fi = 1,2, ... , M, all other being zero. In 
other words, this is just M non-interacting selfdual (critical) Potts models. 

The other special point is 

Xo = X, = 1, (all other xi = 0). (3.18) 

In terms of the original couplings this means K, = 1 + n^=i q}/^^ ^ill other being zero. 
In this case, the M models effectively couple so as to form a single critical Potts model 
with n^li Qfj. states. 



3.2 Two models 



Let us briefly show how to recover the result for M = 2 |]TD[ from the compact formulation 
of Eq. (|3J[5|) . 

Introducing the shorthand notation 5^ = 5((Tj-''\ o"]''"*) for /x = 1,2, the Hamiltonian 
^ reads in this case 



Hij = K,oSi + Ko,S2 + K,,6i62. 



(3.19) 



From Eq. ( p.7| ) we have 

^/qlx,o = b,o = e^'° - 1 

^Xo. = bo. = e^°--l (3.20) 

y^x„ = b„ = e^-'+K'^+Ko. _ (^gX.o ^ gi^o.^ ^ 1^ (3 21) 

and using the selfduality criteria and X,, = Xoo = 1 we readily find the solutions 



1 + fe^--l 



-K.o 



(3.22) 
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3.3 Three models 



For the case M = 3, the Hamiltonian (|3.2| ) reads 



Hi 



K,ooSi + K 0,062 + K 00,63 + 



Ko,,S2S3 + K,o,6i63 + K,,oSi62 + K,,,6i62S3 



(3.23) 



From Eq. ( p.7|) we have 



b 

bo,, 
b,,. 



-.K00m~\~Komm'^Kmo»~\'Kmmo'^K0oo-\-Komo'^Koom 



^g^O««~l~^0«0~l~^00« _|_ g^«o«H~^«oo~h^oo« _|_ g^««o ~l~^«oo ~l~^o«o ^ _|_ 

(e^*°° + e-^"*" + e^°°*) - 1, 



(3.24) 



where the first and second hne each represent three equations that can be obtained by 
cyclically rotating the layer indices. 

Imposing the selfduahty criteria s/qibo,, = y^92?3^»oo and b,,, = ^/qlq2<E we obtain 
after a little algebra 



Mo.. 



^/ g2g3 ^K.oo 1^ + g^o.o _j_ qKoo. 2 



qKo.o -\-Koo. 



-^^ • o o ~l~-^o • o ~l~ -?^o o • 



A 



B 



A 

B ^ 

(e^— - 1) gags + (e^°*° - i) gigs + (e^°°* - i) gig2 + 

(e^— + e^°-° + e^°°* - 2) ^^I^ + gigsgs 
(e^°'° + e^°°* - 1) VgT + {e^''°° - l) ^/q^ x 
(e-^*°° + e^°°' - 1) + (e-^°*° - l) ^/glg^j x 



ig2 



(3.25) 



These expressions generalise those given in Ref. 
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It should be clear that for higher M, manipulating such expressions by direct use of 
the Fourier method becomes extremely cumbersome. 



4 Transfer matrix algorithm 

In Ref. [0 it was shown that the transfer matrix T{L) for coupled Potts models on semi- 
infinite strips of width L may be written in a variety of ways. For a given L, the smaller 
the dimension of T{L) the more efficient will be the computations, since both time and 
memory consumption increase roughly linearly with dim T{L). 

The best choice turns out to be to write the transfer matrix for the loop model on the 



medial lattice Cm IHl, which was referred to as algorithm alg4 in p|. This also gives us 
the advantage of having simple duality relations, cf. Eq. ( p.l5| ), and to be able to treat 
the numbers of states q'2, ijs} as continuous parameters. 



Let us recall from Eqs. ( |3.13|) and (|3.16| ) that the partition function for M = 3 coupled 



models on the self dual manifold can be written, up to a trivial multiplicative constant. 



as 



y _ \^ ^B.oo+Bo,,„Bo,o+B,o,„Boo,+B,,oLi/2 L2/2 L3/2 x 
— 2_^X^ X2 Qi q2 ^3 [^.i-j 

Q 

where we have put xi = x,oo, X2 = Xo,o and X3 = Xoo» for brevity. For simplicity, we 
shall take the lattice C to be the square one, so that £m is once again a square lattice. 

In the loop picture, the symbols • and o refer to the two different ways of splitting 
the vertex, with a definition that alternates between the even and the odd sublattice; 
see Fig. |I|. But thanks to the selfduality, Eq. (|7T) is invariant with respect to a global 



colour conjugation • ^ o, and so one might as well forget about the distinction between 
the sublatticesQ. The number of occurrences of a given vertex splitting (and hence the 
B's) can be counted locally, and are thus easily realised as local Boltzmann factors in 
the transfer matrix. 



^Using the language of Ref. |14|, coupled Potts models constitute a staggered vertex model, which 
however becomes homogeneous exactly on the selfdual manifold. 
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even site 



odd site 



Figure 1: The relation between the random cluster model on the square lattice and 
the loop model on the corresponding medial graph. The clusters consist of connected 
components of coloured edges (thick lines) or isolated sites (filled circles). Loops on the 
medial graph (dashed lines) are defined by the convention that they wrap around the 
cluster boundaries. 

The in ( [4.1| ) count the number of closed loops in each layer fi = 1,2, 3. Despite 
of the non-local nature of the loops, these quantities do not obviate the construction of 
the transfer matrix. Rather, they can be counted locally by writing the transfer matrix 
in the basis of Catalan-like connectivities, as described in 0. Formally, if Ck designates 
the space of pairwise Catalan connectivities among a set of L = 2k points, T{L) acts on 
the product space = Ck^Ck ^Ck- One has dim = k\{k+iy. - 

Weighing the loops with a Boltzmann factor y/g^ that depends on the layer is a 
trivial modification of [Q. An important difference, however, is that the layers are now 
distinguishable. Accordingly, one has simply dim = (dim C^)^. 

As in Ref. we transfer along one of the main directions of the square lattice Cm, 
with periodic boundary conditions in the transverse direction. To ensure that Ck be well 
defined, we must take the strip width to be even, L = 2k. 

We shall also need to consider the situation where all vertex weights x tend to infinity 
in fixed ratios. Defining X2 = ^ and = ^ as the relevant ratios, ( ^.1| ) can be rewritten, 
once again up to an irrelevant multiplicative constant, as 



where the symbol Q' means all colouring figurations in which the local colourings 000 
and • • • do not occur. In Ref. 0], the special case (2:2, 2:3) = (1, 1) with gi = ^2 = Q3 




(4.2) 
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was identified as the non-trivial critical fixed point for three identical coupled models. 



5 Numerical results 

Using sparse matrix factorisation techniques, we have been able to numerically compute 
the first few eigenvalues of the transfer matrices T{L) for Z ( [4. 1| ) and T'{L) for Z' ( [4.2|) for 



even strip widths up to Lmax = 12. The largest matrices had dimension dim T(Lmax) = 
(132)^ but the sparse matrix factorisation necessitates the use of intermediate states 
with L + 2 dangling loop segments, and so the largest sparse matrices involved were of 
dimension (429)'^. 

5.1 Phase diagram and central charge 

As a first test of our algorithm, we checked that it gave the correct eigenvalues at the 
special points {xi,X2,X3) = (1,1,1) and (0,0,0), cf. Eqs. (|3.17| ) and ( p.l8|) . 



We then extracted the (effective) central charge from the leading eigenvalue Aq in the 
standard way [ITBl O: 



/o(^) = /o(oo)-^ + --- (5.1) 

with fo{L) = — log Aq. Using the c-theorem [§], we could readily establish the topology 
of the RG fiows: 

• In the space (xi, 0:2, X3), there is a one-parameter curve along which the partial 
derivatives of the effective central charge with respect to the two perpendicular 
directions vanish identically. Moreover, the corresponding second derivatives are 
strictly negative, so that this curve acts as a "mountain ridge" for Ces- 

• The curve passes through the points (0, 0, 0) and (1, 1, 1), and then goes to infinity 
with some fixed ratios {x2,x'^) = {x2,x'^) that depend on {qi,q2,<l3)- 

• For various values (q'l, (?2, 93) £ [2,4]^, with at most one qa = 2, we observe that 
Ceff is a monotonically decreasing function when going along the curve from the 
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L 



(CefT)max Ac 



2,4,6 0.8523(1) 

4,6,8 0.8476(1) 

6,8,10 0.8481(1) 

8,10,12 0.8483(1) 



0.7133(1) 
0.7136(1) 
0.7076(1) 
0.7061(1) 



1.99348414 



1.93540138 



1.98040147 



1.99047341 



0.00443 



0.00352 



0.00416 



0.00058 



Table 1: Maxima of the effective central charge for the system A^345, and the deviation 
from the value at the decoupling point. 

point (1, 1, 1) towards either (0, 0, 0) or xi(l, x'^), xi oo. In the former case, 
the decrease is rapid and gets more pronounced with increasing system size L, 
signalling the first-order nature (Ceff = 0) of the phase transition in the {qiq2(l3)- 
state Potts model. In the latter case, the decrease is very slight, distinguishing the 
point ( candidate for the non-trivial fixed point of three coupled Potts 

models in the lattice realisation. 

Based on this evidence, we switched to the matrix T'{L) in order to accurately locate 
the point {x2,x'^) and study its critical properties. Invoking again the c-theorem P], 
this was done by scanning the space {x2,x'^) for a maximum in c^s for different system 
sizes. As usually 0, we obtained the best precision by including a non- universal 
term in ( ^.1|) , and so our fits are based on three consecutive strip widths L. The obtained 
maximum can be interpreted as a finite-size pseudo-critical point, which tends to (xg*, x'^) 
as L ^ oo. 

We have concentrated the main part of our computation time on two important special 
cases: A tricritical Ising model and a three-state Potts model coupled with either an Ising 
model or with a tricritical three-state Potts model. In the language of minimal models, 
we refer to these two situations as AI345 = M3 x Ma x M5 and Mi^e = -^4 x M5 x Me 
respectively. 

The corresponding finite-size estimates for {x2,x'^) and for (ccfT)max are given in Ta- 
bles H and||. Rather than directly extrapolating (cefr)max to the L — > 00 limit, we consider 
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L 


X2 


X3 


(Ccff )max 


Ac 


2,4,6 


0.9230(1) 


0.8679(1) 


2.27445933 


-0.0115 


4,6,8 


0.9218(1) 


0.8661(1) 


2.32198052 


-0.0185 


6,8,10 


0.9213(1) 


0.8645(1) 


2.33129236 


-0.0190 


8,10,12 


0.9213(5) 


0.8645(5) 


2.333272 


-0.0225 



Table 2: Maxima of the effective central charge for the system A^456, and the deviation 
from the value at the decoupling point. 

instead for any L its deviation Ac(L) with respect to the corresponding value found at 
the point where the three Potts models decouple. Since this deviation is numerically 
small, it is reasonable to expect that the finite-size corrections to c at {x2-,x'^) are simi- 
lar to those at the decoupling point. Tables |l] and ^ confirm that this is indeed the case: 
the estimates for Ac only depend very weakly on L. As L —>■ 00, we obtain finally the 
following extrapolated values for Ac and for the critical couplings: 

Model M345 ■■ (4% 2:3 ) = (0.8485(1), 0.7053(5)) Ac = -0.0050(3) 
Model M456 ■■ (4% 2:3 ) = (0.9210(5), 0.862(2)) Ac = -0.025(3) 

The values of Ac are seen to compare quite favourably with the two-loop result of the 
perturbative RG (p.24|) , whose numerical values read respectively 



Ac(A^345)rg = -0.0041 Ac(A^456)rg = -0.0175 (5.3) 

Actually, assuming that the RG series is alternating, the fifth-order term not present in 
( p.24| ) would supposedly lead to even better agreement with the numerical results. 

5.2 Higher exponents in the even sector 

Examining the higher eigenvalues of the transfer matrices, one has also access to the 



scaling dimensions by using the formula [ffS 



ML)-ML) = ^-!^^^^ + ... (5.4) 
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We here work in the even sector of the transfer matrix, and besides the identity operator 
(the free energy) we expect to find the various energy-hke primaries. In particular we 
are interested in the primary operators el, discussed in section ^ cf. Eq. (|1.5|) . 

According to Eqs. (CT)-(|2A^) we have < A(e;),A(£^) < 1 and 1 < A{el) < 2 for 



small enough {e^}, and so we expect the first three gaps of the transfer matrix to be 
associated with these operators. 

As in Ref. the scaling dimensions of the operators have been extracted by adding 
to Eq. ( |5.4| ) a non-universal term, and so the fits are based on two consecutive strip 
widths L. 

We consider first the case of coupling three identical models. 

At the non-trivial fixed point {x'2,x'^) = (1,1) the second and third eigenvalues of 
the transfer matrix are degenerate. These eigenvalues are associated with the operators 



^Ai, of ( |1.4|) which generate the two-dimensional irrep of S3, as discussed in the 
Introduction. 

It is important to notice that in the transfer matrix algorithm of the three layers 
were treated symmetrically and only the operators which are symmetric under permuta- 
tion of the layer indices were accessible. In that case, the degenerate eigenvalues asso- 
ciated with the energy doublet were not present in the spectrum of the transfer matrix. 
The fact that we find them here is a strong confirmation of the antisymmetric nature of 
the corresponding operators. 

In Table ^ we report numerical values of the dimension A(£:a) = A(eAi) = ^(^Aa) in 
the case of three coupled Potts models with g = 3. The final result A(£:a) = 0.63(3) is 
in reasonable agreement with the RG value and also agrees well with the Monte Carlo 
result A(£a) = 0.63 ± 0.04 given in 0. 

The fourth eigenvalue of the transfer matrix was already found in and it corre- 
sponds to the symmetric energy operator es of Eq. ( |1.3D . It is a non-trivial check of our 
transfer matrices that the corresponding eigenvalue is identical to that of 0. 

We now turn our attention to the case of three different coupled models. 
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L 


A(eA) 


4,6 


0.6768 


6,8 


0.6648 


8,10 


0.6525 


L — s> oo 


0.63(3) 


RG result 


0.72 



Table 3: Dimension of the antisymmetric energy operator at {x'2ix'^ = (1, 1) ^oi three 
coupled three-state Potts models. 

In the space {x^iX'^) of the couplings we have found that for any given system size 
L there is a unique point (x2,a;3)deg, at which the second and third eigenvalues are 
degenerate. We have determined the location of this point for the models M-zib and 
with L = 4,6, 8, 10 and it is shown on Figs. 0-0. 

In contradistinction to the case of three coupled identical models, the points {x2, X3 ) 
and {x2, a;3)deg do not coincide for finite systems. It is however evident from Figs. 
that these two points become closer as the system size increases. This suggests that in 
the thermodynamic limit they actually coincide. If this scenario holds true, we have at 
the critical point that A{e^) = A{e2). 

Since the finite-size effects in the sequence of points {x'2, a;3)dcg are much stronger than 
is the case for (xg*, Xg*), it is reasonable to measure numerical values of the dimensions at 
the best available estimate for the latter point. We display the corresponding results for 
A{el), A(£:2), A{el) in Tables §-0, once again for the models A^345 and A^456 respectively. 
The numerical results are in quite good agreement with the two-loop RG calculation of 
section ^. Notice that in the case of ^.{62), A.{e^) we compare with the RG values with 
the splitting ignored. 

In Tables we have given the results for ^.{82) and A.{e'^) separately, although we 
do in fact believe that the two dimensions coincide. This point of view is corroborated by 
the fact that the difference between the dimensions of 62 and decreases rapidly with the 
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0.78 
0.77 
0.76 
0.75 
0.74 
0.73 
0.72 
0.71 
0.7 



DEGENERACY 
CENTRAL CHARGE 
EXTRAPOLATION 



0.845 0.85 0.855 0.86 0.865 0.87 0.875 



Figure 2: The location of the critical point for the model M.345, as determined from 
the maximum of the effective central charge. Filled triangles show the positions of the 
maxima {x2, x'^) as found from three-point fits with system sizes ranging from L — 2,4,6 
to L = 8, 10, 12. The asterix indicates the extrapolation of these points to the thermody- 
namic limit, L — > 00. Filled boxes give the positions {x'2, a^3)deg of the degeneracy between 
the second and third eigenvalues of the transfer matrix, with system sizes ranging from 
L = 4 to L = 10. This latter sequence appears to converge to the same hmiting point as 
the central charge data, in support of the symmetry restoration scenario. 
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0.895 I — 

0.89 ■ DEGENERACY 

A CENTRAL CHARGE 

0.885 

M EXTRAPOLATION 

0.88 
0.875 

0.87 
0.865 

0.86 

0.92 0.922 0.924 0.926 0.928 0.93 0.932 0.934 0.936 

Figure 3: Same data as in Fig. ||, but here for the model A^456- 



L 




A(4) 




4,6 


0.75734 


0.7677 


1.2991 


6,8 


0.759 


0.763 


1.2803 


8,10 


0.7564 


0.75647 


1.2761 


RG resuh 


0.746 


0.771 


1.287 



Table 4: Dimensions of the energy operators at {x2,x'^) for the A^345 model. 
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L 








4,6 


0.6849 


0.6857 


1.2958 


6,8 


0.6731 


0.6756 


1.2765 


8,10 


0.6621 


0.6637 


1.2745 


RG result 


0.701 


0.728 


1.43 



Table 5: Dimensions of the energy operators at {x2,x'^) for the M.im model. 

lattice size. For L = 8, 10 the numerical results give Agpiittmg = — ^(^2)1 < 0.0003 

for the AI345 model and Agpntting < 0.002 for the AI456 model. 

In order to be sure that the very small values for Agpiitting obtained at finite size are 
not somehow accidental, we have determined an upper bound for Agputting by directly 
extrapolating the difference between the second and third eigenvalues of the transfer 
matrix at the point {x'^, x'^) . The result is shown below together with the extrapolated 
values for the dimensions of the three energy operators: 



(5.5) 



Model M345 : A(£*) = 1.269 ± 0.002 

A(e*), A(£*) = 0.750 ± 0.005 Asputtmg < 0.003 
Model Ma5q : A(£*) = 1.272 ± 0.002 

A(£*), A(£*) = 0.645 ± 0.003 Asputting < 0.002 

The difference Agputting is well below the one predicted by the RG calculation, which reads 
respectively 0.025 and 0.027 for the models M^Ab and M^m- The numerical work thus 
provides clear evidence that at the non-trivial fixed point the splitting of the dimensions 
of £3 and £2 is actually zero. 



6 Discussion 

In this paper we have shown that coupling M = 3 different Potts models (with ?3 > 
2 and not too large) one obtains a unique non-trivial critical point. The critical properties 

of this point, in particular its central charge and the values of various energetic scaling 
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dimensions, have been determined quite accurately by a perturbative RG treatment and 
found to be consistent with large-scale numerical simulations. An exception is however 
the RG prediction that at two-loop order the degeneracy between the two antisymmetric 
energy operators should be lifted: this prediction has here been discarded on the basis 
of numerical evidence. 

An extension of our investigation to the case of spin-like operators will be published 
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This work forms part of a larger project 0, in which we examine the possible uni- 
versality classes of coupled Potts models, and eventually their relation to the random- 
bond Potts model. In particular, substantial evidence has been accumulated that in 
the random-bond case replica symmetry is not broken, and one can thus hope to make 
analytical progress by studying the unitary models that result from coupling a certain 
number of minimal models. 

We believe that the symmetry properties of the coupled models play an essential role. 
It is in the light of this belief that the present work appears to be interesting: even in 
the absense of an explicit Sm symmetry in the initial action, this symmetry appears 
to be restored at the non-trivial critical point. As far as a putative CFT classification 
of Sm symmetrical critical points goes, we therefore see that the number of models to 
be classified is potentially very large. Indeed, assuming the conclusions of the M = 3 
case to carry over to a general number M of coupled models, one may expect a distinct 
Sm symmetric universality class to arise from coupling any different set of M minimal 
models. If this is true, it would call for a substantial number of new CFTs endowed with 
extended symmetries. Further research along these lines is currently in progress. 
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